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ABSTRACT 

Fewer giants planets are found around M dwarfs than around more massive stars, and this dependence of planetary characteristics 
on the mass of the central star is an important observational diagnostic of planetary formation theories. In part to improve on those 
statistics, we are monitoring the radial velocities of nearby M dwarfs with the HARPS spectrograph on the ESO 3.6 m telescope. We 
present here the detection of giant planets around two nearby M0 dwarfs: planets, with minimum masses of respectively 5 Jupiter 
masses and 1 Saturn mass, orbit around Gl 676A and HIP 12961. The latter is, by over a factor of two, the most massive planet 
found by radial velocity monitoring of an M dwarf, but its being found around an early M-dwarf is in approximate line with the 
upper envelope of the planetary vs stellar mass diagram. HIP 12961 ([Fe/H]=-0.07) is slightly more metal-rich than the average solar 
neighborhood ([Fe/H]=-0.17), and Gl 676A ([Fe/H=0.18) significantly so. The two stars together therefore reinforce the growing 
trend for giant planets being more frequent around more metal-rich M dwarfs, and the 5 Jupiter mass Gl 676Ab being found around 
a metal-rich star is consistent with the expectation that the most massive planets preferentially form in disks with large condensate 
masses. 

Key words. Stars: individual: Gl 676A - Stars: individual: HIP 12961 - Stars: planetary systems - Stars: late-type - Techniques: 
radial-velocity 



1. Introduction 

Much recent theoretical work has gone into examining how 
planet formation depends on stellar mass, because stellar mass 
significantly changes the physical conditions which control the 
formation of planets. A comparison, for instance, of the planet 
populations around Sun-like stars on one hand, and around M 
dwarfs on the other hand, probes the sensitivity of the plane- 
tary formation process to several physical parameters: around 
lower mass stars gravity (hence disk rotation speed), tempera- 
ture (which regulates the position of the ice line) are both lower, 
and, perhaps most importantly, disk mass scales approximately 
linearly with stellar mass (e.g. Scholz et al.|2006). 



Within the "core acc retion" paradigm, |Laughlin et al. ( 2004| l, 
Ida & Lin| pOOS] ), and |Kennedy & Kenyon| ( |2008) > all predict 
that giant planet formation is inhibited around very-low-mass 
stars, while Neptune-mass planets should inversely be com- 
mon. Within the same paradigm, but assuming that the prop- 
erties of protoplanetary disks, contrary to observations, do not 
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change with stellar mass, |Kornet et al.| ([2006 ) predict instead 
that Jupiter-mass planets become more frequent in inverse pro- 
portion to the stellar mass. Finally, Boss] (2006) examines how 
planet formation depends on stellar mass for planets formed by 
disk instability, and concludes that the frequency of Jupiter-mass 
planet is largely independent of stellar mass, as long as disks are 
massive enough to become unstable. One needs to note, though, 
that proto-planetary disks of a realistic mass are likely to be grav- 
itationally stable out to at least 10 AU. Planets can thus form 
through gravitational instability only beyond that distance, in 
a separation range only skimmed by radial velocity monitoring 
and probed mostly by microlensing searches and direct imaging. 
Massive planets formed by gravitational instability and found 
well within 5 AU must thus then have migrated inward. How 
giant planets migrating in the massive disks needed for gravita- 
tional instability can escape accreting enough mass to become a 
brow n dwarf (> 13 M Jup ) is unclea r Stamatellos & Whitworfh 
(e.g. [2009]); [Kratter et al.|(e.g.|20T0|> 



Observationally, just a dozen of the close to 400 planetary 
systems currently known from radial velocity monitoring, are 
centered around M dwarfs (M<0.6 M Q ) Q This no doubt re- 
flects in part a selection bias, since many more of the intrisically 
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Parameter 

Spectral Type 

V 

J 

K s 

71 

Distance 

My 

Mk„ 
M b „i 
L* 

v sin i 
[Fe/H] 



[mas] 
[pc] 



[L s ] 
[km s~ 



Gl 676A 

MOV 
9.585 ± 0.006 
6.711 ±0.020 
5.825 ± 0.029 

2.73 
60.79 ± 1.62 
16.45 ± 0.44 
8.50 ± 0.06 
4.74 ± 0.06 
7.47 
0.082 
1.6 ± 1.0 
0.18 
0.71 



HIP 12961 

MOV 
10.31 ±0.04 
7.558 ±0.021 
6.736 ±0.018 

2.61 
43.45 ± 1.72 
23.01 ±0.91 
8.50 ± 0.09 
4.93 ± 0.09 
7.54 
0.076 
1.5 ± 1.0 
-0.07 
0.67 



Notes 

~W~ 
(b) 
(c) 
(c) 
(d) 
(c) 



(f) 
(g) 



Table 1. Observed and inferred stellar parameters for Gl 676A 
and HIP 129 61. Notes: (a) |Hawley etalT1([T996)l fo r Gl 676A, 
and |Luyten| ( fl980] l for HIP 12961; (b) |Koen et al.| ( |2002> for 
Gl 676A, and co mputed from the TYCH O (B r ,V r ) photometry 
for HIP 12961 ; (c)|Skrutskie et al.|([2006); (d) computed from J- 



K s with th e Leggett et al. ^OOTjBC^ vs J-K s relation; (e) 
Leeuwen (2007); (f) computed from Mk and V - K s using 



van 

"Ee 



Schlaufman & Laughlin (2010) calibration; (g) computed from 



M^ s using the Delfosse et al.| (2000) calibration; both masses are 
at the upper edge of the validity range of that calibration, and 
they might therefore have somewhat larger errors than its 10% 
dispersion. 



brighter solar-type stars than of the fainter M dwarfs have been 
searc hed for planets, but there is increasing statistical evidence 
Bonfils et al-pOoS} |Endl et al.||2006l |Johnson et alj|2"007 



(e.g. 



2010) that M dwarfs also genuinely have fewer massive planets 
(~Mj up ) than the more massive solar-type stars. They may, on 
the other hand, and though no rigorous statistical analysis has yet 
been performed for that planet population, have a larger preva- 
lence of the harder to detect Neptune-mass and super-Earth plan- 
ets: a quarter of the ~30 planets with M sin(i) < 0.1My„ ; ,known 
to date orbit an M dwarf, when solar-type stars outnumber 
M dwarfs by an order of magnitude in planet-search samples. 
Conversely, the highest mass planets known around M dwarfs 
are the M sin(i) = 2M y „ p Gl 876b (|Delfosse et al.|[T998| |Marcy 



|etal.|1998| ) and HIP7943 lb ( |Apps et al.|2010| ), and atalarger or- 
bital separation of ~3 AUs t he M ^^3.5Mj I , p OG LE-2005-BLG- 
071 Lb microlensing planet (Dong et al. 2009), when over two 
dozen planets with masses over 10My Hp are known around solar 
type stars. The statistical significance of that difference however 
remains modest, since the M dwarfs searched for planets only 
number in the few hundreds, when the apparent fraction of these 
very massive planets is under 1% around solar type stars. 

We present here the detection of two giant planets around 
M0 dwarfs, a M sin(i) = 0.354M 7t(p planet around HIP 12961, 
and a M sin(i) = 4.87M y „ p planet around Gl 676A. 

2. Stellar characteristics 

Table [TJ summarizes the properties of the two host stars, which 
we briefly discuss below. 

2.1. HIP 12961 

HIP 12961 (also CD-23degl056, LTT 1349, NLTT 8966, 
SAO 168043) was not identified as a member of the 25 pc vol- 
ume until the publication of the |Perryman & ESA] ( |1997| ) cat- 



alog, and has attracted very little attention: it is mentioned in 
just 4 literature references, and always as part of a large catalog. 
Because HIP 1291 does not figure in the Gliese & JahreiB (1991) 
catalog, it was omitted from the Hawle yet al.| ( |1996[ ) spectral at- 
las of the late-type nearby stars. SIMBAD shows an M0 spectral 
type, which seems to trace back to a classification of untrace- 
able pedigree listed in the NLTT catalog ( |Luyten|fl980| ), while 
Step henson] (fT986 ) estimated a K5 type from low-dispersion ob- 
jective prism photographic plates. The absolute magnitude and 
color of HIP 12961, M y =8.50 and V-K=3.57, suggest that it s 
older NLTT spectral type is closer to truth (e.g. Leg gett|1992[ ). 
We adopt this spectral type for the reminder of the paper, but note 
that a modern classification from a digital low resolution spec- 
trum is desirable. HIP 12961 has fairly strong chromospheric 
activity, with 90% of stars with spectral types K7 to Ml in 
the HARPS radial velocity sample (which however reject the 
most active stars) having weaker Can H and K lines, and just 
10% stronger lines. The 2MASS photometry (Table [T]t and the 



Leggett et al. (2001 1 J - K colour vs bolometric relation result 
in a K-band bolometric correction of BCk = 2.61, and together 
with the parallax in a 0.076 L luminosity. 



2.2. Gl 676A 

The Gl 676 system (also CCDM J17302-5138) has been rec- 
ognized as a member of the immediate solar neighborhood for 
much longer, figuring in the original Gliese (1969) catalog of 
the 20 pc volume. It consequently has 15 references listed in 
SIMBAD, though none of those dedicates more than a few sen- 
tences to Gl 676. The system comprises Gl 676A (also CD- 
51 10924, HIP 85647, CPD-51 10396) an d Gl 676B, with re- 
spective spectral types of MOV and M3V (Ha wley et al.|[l 996) 
and separated by ~50" on the sky. At the distance of the system 
this angular distance translates into an ~800 AU projected sepa- 
ration, which is probably far enough that Gl 676B didn't strongly 
influence the formation of the planetary system of Gl 676A. 
Gl 676A is a moderately active star, with a Can H and K emis- 
sion strength at the third quartile of the cumulative distribution 
for stars with spectral types between K7 and Ml in the HARPS 
radial velocity sample. The 2MASS photometry (Table [TJ and 
the |Leggett et al.| ( [2001| ) J - K colour-bolometric relation result 
in a K-band bolometric correction of BCk = 2.73, and together 
with the parallax in a 0.082 L luminosity. The |Delfosse et aLj 
(2000) K-band Mass-Luminosity relation results in masses of 
respectively 0.71 and O.29M for Gl 676A and Gl 676B. The 
former is at the edge of the validity range of the Delfosse et al. 
(2000) calibration, and might therefore have somewhat larger 
uncertainties than the ~10% dispersion in that calibration. 



3. HARPS Doppler measurements and orbital 
analysis 

We obtained measurements of Gl 676A and HIP 12961 with 
HARPS (High Accuracy Radial velocity Planet Searcher Mayor 



eTaT1[2003 



as part of the guaranteed-time program of the in 
strument consortium. HARPS is a high-resolution (R = 115 
000) fiber-fed echelle spectrograph, optimized for planet search 
programs and asteroseismology. It is the most precise spectro- 
velocimeter to date, with a long-term instrumental RV accuracy 
well under lms" 1 



(e.g. 



Lovis et al. 2006 



Mayor et al. 



2009) 



When it aims for ultimate radial velocity precision, HARPS uses 
simultaneous exposures of a thorium lamp through a calibration 
fiber. For the present observations however, we relied instead 
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Fig. 1. HARPS radial velocities of HIP 12961 as a function of 
barycentric Julian day (top panel), spectral power window (mid- 
dle panel) and Lomb-Scargle periodogram of these velocities 
(bottom panel). The horizontal lines mark false alarm proba- 
bilities equivalent to 1, 2 and 3 sigmas significance levels for 
Gaussian noise. 



on its excelent instrumental stability (nightly instrumental drifts 
< 1 m s ). Both targets are too faint for us to reach the stability 
limit of HARPS within realistic integration times, and dispens- 
ing with the simultaneous thorium light produces cleaner stellar 
spectra, more easily amenable to quantitative spectroscopic anal- 
ysis. The two stars were observed as part of the volvme-limited 
HARPS sea rch for planets (e.g. |Moutou et~aL]|2009[ |Lo Curto 
et al.||2010| . While generally refered to as F-G-K stars, for the 
sake of concision, the targets of that program actually include 
M0 dwarfs ( |Lo Curto et al.|2010) . 

We used 15 mn exposures for both stars, obtaining median 
S/N ratios (per pixel at 550 nm) of 53 for the V=9.58 Gl 676A, 
and 49 for the V=10.31 HIP 12961. The 69 and 46 radial veloc- 
ities of Gl 676A and HIP 12961 (Tables [3] and g] only available 
electronically) were obtained with the standard HARPS reduc- 
tion pipeline, based on cross-correlation with a stellar mask and 
on a precise nightly wavelength calibration from ThAr spectra 
( |Lovis & Pepe|200"7] l. The median internal errors of these veloc- 
ities are respectively 1.9 and 2.8 ms _1 , and include a ~ 0.2 ms" 1 
noise on the nightly zero-point measurement, a ~ 0.3 ms _l un- 
certainty on the instrumental drift, and the photon noise com- 
puted from the full Doppler information content of the spectra 
( Bouchy et al.|200T i. The photon noise contribution completely 
dominates the error budget for these moderately faint sources. 
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Fig. 2. HARPS radial velocities of HIP 12961 phased to the 
57.44 days period, overlaid with the adjusted Keplerian orbit 
(top panel). Residuals from the Keplerian orbit as a function of 
barycentric Julian day (middle panel), and Lomb-Scargle peri- 
odogram of these residuals (bottom panel). The horizontal lines 
mark false alarm probabilities equivalent to 1, 2 and 3 sigmas 
significance levels for Gaussian noise. 



3.1. A Saturn-mass planet around HIP 12961 

The computed radial velocities of HIP 12961 vary with a 
~ 60m s~ l peak to peak amplitude (Fig.[T| top panel), an order of 
magnitude above their 2.6 m s average photon noise, and well 
above the <10 ms~ l maximum jitter expected from the chro- 
mospheric activity. The variations show no correlation with the 
bisector span (rms 6m s ), the depth (10.65%, with 0.07% rms) 
or width (3.566 km s~ l , with 14m s rms) of the correlation pro- 
file, or any of the standard stellar activity diagnostic, making or- 
bital motion by far their most likely cause. The Lomb-Scargle 
periodogram of the velocities shows one highly significant peak 
at 57.45 days (Fig.[T[ middle and bottom panel), as well as its 4 
aliases at +-1 sidereal and civil days. Phasing of the velocities 
on that period shows well sampled smooth variations (Fig.[2]top 
panel). A Keplerian fit (Table |2]i yields a moderately eccentric 
orbit (e = 0.2) with a 25 ms -1 semi-amplitude. 

The rms amplitude of the residuals from that orbit is 
3.8 m s _1 (Fig. [2] middle panel), significantly above the 2.6 m s _1 
average photon noise. The square root ot the reduced^- 2 of the fit 
is consequently 1.5. The radial velocities may therefore contain 
information beyond the detected planet, but the highest peak in 
the periodogram of the radial velocity residuals (16.6 days) only 
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has 2 cr significance. It also coincides with a signal in the pe- 
riodogram of the correlation profile's depth, and it is broadly 
consistent with the stellar rotation period expected from the sig- 
nificant chromospheric activity. This peak, if not just noise, is 
therefore much more likely to reflect rotational modulation of 
stellar spots than a planet. There is no current evidence for addi- 
tional planets in the system. 

Together with the 0.67 M (Table [TJ stellar mass, the or- 
bital parameters translate into a minimum companion mass of 
0.35Mj up , or 1.2 Saturn-mass, with a 0.25 AU semi-major axis. 
For an Earth-like albedo of 0.35 the equilibrium temperature at 
that distance from a 0.076 L Q (Table[T]i luminosity star is 263 K, 
slightly higher than the terrestrial 255 K but below the ~270 K 
threshold for triggering a runaway greenhouse effect |Selsis et al. 
(2007 ). A putative massive moon of HIP 12961b could therefore 
be potentially hospitable to life. 

The a priori geometric probability that HIP 12961b transits 
across HIP 12961 is approximately 1%. As usual for planets of 
M dwarfs, the transit would be deep (~2.5%) and therefore well 
suited to high quality transmission spectroscopy of the plane- 
tary atmosphere, as well as to searches for transits by plane- 
tary moons. This high potential return offsets the long odds to 
some extent, and the deep transits would be within easy reach 
of amateur-grade equipment. The star will thus be well worth 
searching for transits, once additional radial velocity measure- 
ments will have narrowed down the time windows for potential 
planetary transits. 

3.2. A massive long period planet around Gl 676A 

The computed velocities of Gl 676A (Table|3]l exhibit unambigu- 
ous variations of several hundred ms~' with a period slightly 
over our current observing span, superimposed upon a slower 
drift (Fig|3| upper panel). The correlation profile depth (13.61%, 
with 0.12% rms), its width, the bisector span (5 ms~' rms), and 
the chromospheric indices show no systematic variations that 
would correlate with the radial velocity changes. 

A fit of a Keplerian plus a constant acceleration to the ra- 
dial velocities (Table|2]i yields a period of 1056.8 + 2.8 days, a 
semi-amplitude of 122.8 ± 1.9 ms _1 , and a 10.7 ms^'yr -1 ac- 
celeration. 

Together with the 0.71 M Q (Table[TJ stellar mass, the orbital 
parameters imply a companion mass of 4.9 M Jup and a 1.82 AU 
semi-major axis. At the 16.5 pc distance of the Gl 676 system, 
the minimum astrometric wobble of Gl 676A, for a sin(i) = 1 
edge-on orbit, is +0.67 mas. This is within reach of both the 



200 



FGS instrument on HST (e.g. Martioli et al.|2010|) and imagers 


on 8m-class telescopes (|Lazorenko et al. 


2009). The companion 


of Gl 676A therefore belongs to the smal 


group of non-eclipsing 



planets for which the inclination ambiguity can potentially be 
lifted with existing instruments. We have started such astromet- 
ric observations using the FORS2 imager of the ESO VLT. 

The line of sight acceleration of Gl 676A by its 
0.3M G G1 676B stellar companion is of order GMj//^ s , and 
from the 800 AU projected separation is therefore under 
0.1 ms _1 yr _1 . The two orders of magnitude discrepancy be- 
tween the observed acceleration and that expected from Gl 676B 
demonstrates that the system contains an additional massive 
body. That body could be planetary if its separation is under 
-15 AU (0.9"), stellar if that separation is above -40 AU (2.4"), 
or a brown dwarf for intermediate separations. Adaptive optics 
imaging could easily narrow down these possibilities, as will the 
continuing radial velocity monitoring to constrain the curvature 
of the radial velocity trend, and the on-going astrometric effort. 
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Fig. 3. HARPS radial velocities of Gl 676A as a function of 
barycentric Julian day (top panel), spectral power window (mid- 
dle panel) and Lomb-Scargle periodogram of these velocities 
(bottom panel). The horizontal lines mark false alarm proba- 
bilities equivalent to 1, 2 and 3 sigmas significance levels for 
Gaussian noise. 



The rms amplitude of the residuals around the 
Keplerian+acceleration orbit is 3.4 ms~', significantly above 
our 1.7 ms~' average measurement error. The square root of 
the reduced x 2 of the fit is consequently 2.0, indicating that 
the residuals contain structure above the photon noise. The 
highest peak in a Lomb-Scargle ( |Lomb 1976 Scargle 1982) 
periodogram of the residuals however only rises to a level 
equivalent to a 1 cr detection (FigQ. There is therefore no 
immediate evidence for additional planets in the system. The 
excess residuals may simply reflect jitter from the moderate 
stellar activity of Gl 676A, or alternatively they could be 
early signs of multiple additional planets, which additional 
observations would then eventually disentangle. 



4. Discussion 

As discussed above, HIP 12961 and Gl 676A are orbited by 
giant planets with minimum masses of approximately 0.5 and 
5 Jupiter masses. The latter is twice the M sin(i) = 2My„ p of 



Gl 876b pelfosse et al.|1998l|Marcy et al.|1998| l and HIP7943 lb 
(Apps et al. 2 010) 1, previously the highest mass planets found 
by radial velocity monitoring of M dwarfs, and above the 3.8 
or 3.4 My u/ ,(from two degenerate solutions) of the OGLE-2005- 
BLG-071Lb ( [Dong et al.||2009] > microlensing planet. The MOV 
Gl 676A however is significantly more massive (0.71 M Q , 
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Fig. 4. Top panel: HARPS radial velocities of Gl 676A as a func- 
tion of orbital phase, after subtraction of adjusted linear drift and 
overlaid with the adjusted Keplerian orbit. Middle panel: resid- 
ual of the HARPS radial from the adjusted orbit and linear drift, 
as a function of time. Bottom Panel: Lomb-Scargle periodogram 
of the residuals The horizontal lines mark false alarm probabili- 
ties equivalent to 1 and 2 sigmas significance levels for Gaussian 
noise. 

Table 2. Orbital elements for the Keplerian orbital models of 
HIP 12961 and Gl 676 A. 



Table [T} than the M4V Gl 876 ( 0.33 M |Correia et al.| [2010) 
andtheM3VHIP79431 (0.49 M B |Apps et al.|2010) . The higher 
mass of its planet therefore remains in approximate line with the 
current upper envelope of the planetary versus stellar mass dia- 
gram. These most massive planets are rare at any stellar mass, 
with an occurence rate under 1 %, suggesting that they can form 
only under the most favorable conditions. They have been sug- 
gested to form through gravitational instability, with their lower 
mass counterparts forming by core accretion. Proto-planetary 
disks of any realistic mass, however, are expected be gravita- 
tionally stable out to beyond 10 AU. If Gl 676 Ab formed through 
gravitational instability, it would therefore have undergone much 
inward migration, through a very massive disk. How it could es- 
cape accreting enough mass during this migration to become a 
brown dwarf is unclear. 

Gl 676 A and HIP 12961 increase the sample of M dwarfs 
with giant planets (Saturn-mass and above) from 7 to 9, and 
therefore offer an opportunity to evaluate the trend ( |Johnson"& 
Apps 2009 ; Schlauf man & Laughlin|20T0] l for giant planets be- 
ing more common around more metal-rich M dwarfs. Adopting 
the very recent Schlaufman & Laughlin|( |2010[ ) metallicity cali- 
bration of the M^ s vs V-Kj plane, which finds metallicities ap- 



55500 proximately half-way between those of the earlier |Bonfils et al. 



( 2005 1 and Johnson & Apps ( 2009 ) calibrations, the metallicities 
ofGl 676A and HIP 12961 are 0.18 and -0.07. Both values are 
above the [Fe/H] = -0.17 average metallicity for the solar neigh- 
borhood in the Schlaufman & Laughlin|( |2010| l metallicity scale, 
the latter very significantly so. The two new planets therefore 
clearly reinforce the incipient trend, and help suggest that more 
massive planets are found around more metal-rich M-dwarfs. 
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Element 


HIP 12961 


Gl 676A 


y [kms -1 ] 


33.0463 ± 0.0014 


-39.108 ±0.091 


dy/dt [ms~' yr~'] 




10.66 ±0.61 


Epoch [BJD] 




2450000 


P [days] 


57.435 ± 0.042 


1056.8 ±2.8 


e 


0.166 ±0.034 


0.326 ± 0.009 


oj [deg.] 


272 ± 13 


85.7 ± 1.4 


TO [BJD] 


2454428.4 ± 2.0 


2455411.9 ± 3.0 


Kl [ms" 1 ] 


24.71 ±0.86 


129.3 ± 1.2 


a [AU] 


0.25 


1.82 


M sin(i) [M Jup ] 


0.35 


4.9 


cr(O-C) [ms- 1 ] 


3.9 


3.4 


N 


46 


69 


N 

1 * par 
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Table 3. Radial-velocity measurements and error bars for 
Gl 676A All values are relative to the solar system barycenter, 
and corrected from the small perspective acceleration using the 
Hipparcos parallax and proper motion. Only available electroni- 
cally. 
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Table 4. Radial-velocity measurements and error bars for 
HIP 12961 All values are relative to the solar system barycen- 
ter, and corrected from the small perspective acceleration using 
the Hipparcos parallax and proper motion. Only available elec- 
tronically. 
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